Linear Response in Complex Systems: CTRW and the Fractional Fokker-Planck 
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We consider the linear response of systems modelled by continuous-time random walks (CTRW) 
and by fractional Fokker-Planck equations under the influence of time-dependent external flelds. 
We calculate the corresponding response functions explicitely. The CTRW curve exhibits aging, i. e. 
it is not translationally invariant in the time-domain. This is different from what happens under 
fractional Fokker-Planck conditions. 
PACS numbers: 05.40.-a, 02.05.-r 



I. INTRODUCTION 

Many systems, such as polymer chains and networks, 
proteins, glasses and charge-carriers in semiconductors 
are characterized by extremely slow relaxation processes. 
The microscopic mechanisms leading to such slow re- 
laxations differ considerably from system to system, 
but their microscopic manifestations often correspond to 
power-laws. A recently introduced approach to slow re- 
laxation in time-independent external fields is based on 
fractional Fokker-Planck equations (FFPEs) ||l|,p| or on 
fractional Master equations JS). The FFPE describing 
subdiffusive behavior in an external field reads: 



d_ 
dt 



P{x,t) 



qDI ''CppPix^t), 



(1) 



where P(x, t) is the pdf to find a particle (walker) at point 
X at time t. fn Eq.(yj) oDf. '^ is the fractional Ricmann- 
Liouville operator (0 < 7 < 1) defined through ||l| 



,D]-'<Z{t) = 



dt'- 
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Here t = can be associated with the time at which the 
system was prepared. Furthermore, in Eq.(|l|) Cpp is the 
Fokker-Planck operator. 



Cfp = KA- ^Vf (.x) 



(3) 



fn what follows we concentrate on the case in which the 
acting force f{x) = — VC/ is homogeneous and of magni- 
tude E. fn this case we have Cpp — KA — /iEV. 

The fractional Fokker-Planck equation, Eq.(^, has 
turned out to be useful in describing a broad range of phe- 
nomena connected with anomalous diffusion |l|]. There 
are essentially two main mechanisms leading to long-time 
memory in the behavior of complex systems. One is re- 
lated to the hierarchical structure of the modes of the 
system, as is the case for polymer chains and networks Q] 
and for rough interfaces pj . The other is associated with 
the diffusion in complex, (almost) quenched potential 



landscapes, which occur in glassy systems (ranging from 
normal glasses to proteins) . A phenomenological descrip- 
tion of free relaxation (i.e. the time-evolution of a system 
which is prepared in a nonequilibrium initial condition 
and then evolves under time-independent external condi- 
tions) in terms of fractional Fokker-Planck equations is 
reasonable in both situations. Thus, the FFPE, Eq.(|^) 
can be viewed as a phenomenological linear-response the- 
ory for a system with long memory in contact with a 
heat bath ^. This equation can be derived systemat- 
ically from the continuous-time random walk (CTRW) 
scheme using the standard Kramers-Moyal procedure. 
The equivalence of the two approaches was discussed in 
|7| using the subordination property of CTRWs; see also 
|3,||. Thus, the free relaxation properties of a CTRW 
system are closely reproduced by Eq.(Q) and can be ex- 
pressed in terms of Mittag-Leffler functions ||l|,^ . 

In what follows we concentrate on the response of the 
two models to time- dependent external fields, and discuss 
the mean current (velocity) and mean polarization (co- 
ordinate) as a function of time in systems under pulsed 
or sinusoidal external fields. We show that the linear 
response to time-dependent fields predicted by FFPE 
(which is exactly the same as for hierarchical models) 
differs strongly from what is expected under CTRW con- 
ditions. Thus, the response to time-dependent fields of 
the systems described by FFPE is mainly infiuenced by 
the values of the field at times immediately preceding the 
observation time t, while the memory of the earlier his- 
tory fades out. On the other hand, in CTRW-systems 
the response of the current to external fields is local in 
time, while the corresponding susceptibility decays. In 
such systems the polarization depends much on the early 
history of the system. 



II. LINEAR RESPONSE WITHIN THE FFPE 
SCHEME 

Let us consider a system whose dynamics is described 
by a fractional Fokker-Planck equation and concentrate 



on the mean particles' displacement under the action of 
an external force. This mean displacement X(t) is given 
by 



/oo 
xP(x, t)dx. 
-OO 



(4) 



Multiplying Eq. (||) by x and integrating it over the whole 
space we obtain: 



— /xP(x,t)dx= oDI^'^ K /xAP(x,i)dx 
- ^lE{t) I xVP(x,t)dx 



(5) 



The left hand side of Eq.(|) is nothing but ^X(i), 
whereas the right hand side can be simplified by inte- 
gration by parts. Since P(x, t) and its derivatives with 
respect to the coordinates vanish at infinity, the first in- 
tegral vanishes and the second one is unity. Hence 



dt 



X(t) 



oDl-'^^iE{t)■ 



(6) 



Note that such kind of response is typical for complex 
and hierarchically built systems, like polymer chains and 
networks, see Q. 

As a simple example let us consider a chain of iV 3> 1 
beads connected by harmonic springs and immersed in a 
viscous fluid (a Rouse-chain). We assume the first bead 
of the chain to be tagged (say, charged) and experience 
the external field E. The motion of the beads is governed 
by the equation 






xi)+E(t) + fo(t) 



for the first bead (bearing number 0), 



c- 



dt 
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for the last one and 



dt 



(7) 



(8) 
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for aU other beads, < i < iV. In Eqs.(0) to (|) f,(i) are 
Gaussian, (5-correlated forces with zero mean. Note that 
due to the linearity of Eqs.(0) to (Q) the mean positions 
and velocities of the beads (averaged over the realizations 
of the {i{t)) follow equations similar in form, but where 
the random forces fi(t) are omitted. Let us suppose that 
the velocity response of the first bead of the chain (av- 
eraged over the realizations of the random forces fi(i)) is 
described by the memory function M{t). Supposing that 
the initial mean velocities of the beads vanish, this leads 
after a Laplace transformation to 



vo(A) = Axo(A) = Af(A)E(A). 



(10) 



If N is very large, N -^ oo, the chain can be considered 
as infinite; then the subchain starting with bead 1 has the 
same linear response properties as the entire chain start- 
ing with bead 0. Thus, M{t) can be found using Eq.(|j) 
and noting that the equation of motion for bead 1 has 
the same form, Eq.(|lO|), with the acting force now being 
F = fc(xo — Xi) instead of E. The Laplace-representation 
of the equation of motion for the first bead reads: 



Cvo(A) = -k 



vo(A) vi(A) 



A 



A 



E(A) 



(11) 



(where we set xo(A) ~ vo(A)/A and xi — vi(A)/A). 
Moreover, the analog of Eq.(|l(]|) for the velocity of the 
bead 1 reads: 



vi(A) = A/(A)fc 



vo(A) vi(A) 



A 



A 



(12) 



The solution of Eqs.(|lol) to (|l|) gives M{X) = -A/2fc + 
VA2/4fc2 + X/Ck. Thus, for A ^ 0, Af (A) ~ Tq/^VA with 
To = l/C^j which is exactly the Laplace-representation 



-1/2 



of a semi-derivative, Tq'^ qD^^^ with 7 = 1/2 in Eq.(ph. 
This is the value which we will use in our numerical ex- 
amples in what follows. 

In the case where the friction coefficients d and/or the 
spring constants ki differ from site to site, other values 
of 7 may be obtained Q| . The same holds also for more 
complex (higher-dimensional, fractal or tree-like) struc- 
tures. For example, the case when a tagged monomer is 
attached to a membrane corresponds to 7 = 2/3 [ [lO[ . 



III. LINEAR RESPONSE OF A CTRW SYSTEM 

Let us now turn to the linear response in the frame- 
work of the CTRW model introduced by MontroU and 
Weiss ||ll|]. This model was extremely successful in the 
explanation of dispersive transport in amorphous semi- 
conductors |1^; see [ p!3|Jl4| for reviews. As is usual in 
CTRW, we envisage an ensemble of noninteracting par- 
ticles which may be influenced by external fields (say, the 
particles are charged). The particles follow then CTRWs, 
i.e. sequences of jumps. The time intervals ti between 
the jumps are uncorrelated. Of interest are waiting times 
which follow power-law distributions, i.e. 



,p{t) = j/ii+t/Toy+\ 



with < 7 < 1. (13) 



The physical motivation for such '(/'(i)-forms may be ra- 
tionalized usin g ra ndom traps whose energy distribution 
is exponential |15|. In what follows we put tq = 1 and 
work in dimensionless time units. 

A basic quantity in the CTRW formalism is Xn {t) , the 
probability to make exactly n steps up to time t. In the 
standard, decoupled CTRW picture (in which the spatial 
transition probabilities between the lattice sites are in- 
dependent of the waiting-times) under time-independent 



field the probability distribution P(x, t) of finding a par- 
ticle at r at time t given that is started at at time 
obeys |16 



P(x,i)= VP„(x)x„(i), 



(14) 



where Pn(x) is the probability to reach x from in n 
steps. We note that Eq.(|l|) states that the CTRW is a 
random process subordinated to simple random walks un- 
der the operational time given by the Xn(i)-distribution. 
Note that exactly this property is the starting point for 
the derivation of FFPE in Ref. m. Let us now turn to 
the case of time- dependent fields. In this case the simple 
subordination relation, Eq.([l^) breaks down, since Pn(x) 
starts to depend on the actual value of external field at 
time instants of steps, and thus on the actual times of 
steps, and not only on their number. On the other hand, 
the mean velocity or mean displacement of the particles 
can still be easily found. 

Let us discuss the linear response of an ensemble of 
random walkers performing CTRWs to a changing exter- 
nal field. We consider some physically short time interval 
dt. Let dN be the mean number of steps performed dur- 
ing dt. The mean displacement during dt is dX ~ xdN 
where x is the mean displacement per step depending on 
the actual value of the external field E: 



^E = ^x,, 



(x, ■ E) 
ksT 



(15) 



According to the theory of CTRW, x„ (A) , the Laplace- 
transform of Xn{t), reads XnW — V'(-^)" [1 — i^W] /^: 
Ref. jlj]. From Eq.([l8|) we now have that the Laplace- 
transform of /(t) reads: 



/(A) 



V^(A) 
l-^(A) 



(20) 



The Laplace transform of ?A(i), Eq.(p^), for small A is 
known to be V(A) = 1 - A'T(1 -7) O- The in- 
verse Laplace-transform of Eq.(|2^) (for longer times) thus 
reads: 



m = ^^F- 



(21) 



Note that the convergence to the behavior given by 
Eq.(pl|) can be very fast. Moreover, effective interpo- 
lating forms valid both at short and at longer times can 
be obtained. As an example let us consider the case 
when ip{t) is a one-sided (extreme) Levy-law, ip{t) — 
L(i,7,— 7) (with < 7 < 1), whose Laplace-transform 
is a stretched exponential ip{u) = exp(— A'''). This leads 
to /(A) = [exp(A''') — 1] . Note that the behavior for 
large A, i.e. short t, corresponds to that of ip{t), since 
the denominator of /(A) is dominated by the exponen- 
tial term. For small A the asymptotic behavior given by 
Eq.(pl|) sets in. The transition between the two types of 
behaviors takes place at i ~ 1, and thus at longer times 
Eq.(|l]) holds. 

Thus, for t > I one has 



Here ^ is the mobility tensor. The sum in the second 
expression runs over all nearest neighbors vectors, fcs is 
the Boltzmann constant and T is the temperature. We 
obtain now 



dX = fiE{t)dN, 



(16) 



where dN^N{t + dt)-N{t) ^ dtJ2Zo^ft^n{t). Thus, 
the typical particles' velocity (which is proportional to 
the particles' current) is given by: 



dX 



where 



V(t) = ^ = /(i)ME(i). 



/w = E"v^^»w 



i=0 



(17) 



(18) 



Now the current density is J(i) = neV(i), where n is the 
density of charge carriers (assumed to be homogeneous) 
and e is their charge; the analogue of Eq. ([l7|) holds also 
for the current. 



J(i) = /(t)aE(i), 
where a — ne/i is the conductivity tensor. 



(19) 



V(t) 



Sin7r7 



i'^-VE(t). 



(22) 



A similar equation holds, of course, for the current J(i) 
flowing through the system: 



, , , sm 7r7 „ 1 

J(i) = ^F"VE(t). 



(23) 



The current response of a CTRW-system to an exter- 
nal field is local in time and depends explicitly on the 
time elapsed after the system was prepared. Systems in 
which the response to an external agent depends explic- 
itly on the delay between preparation time and measure- 
ment time are referred to as aging systems. This kind 
of behavior is found to be pronounced in CTRWs with 

0<7<lllH§]. 

To obtain the particle's position and therefore the po- 
larization of the medium we simply have to integrate 
Eq.(p2[) over time. Hence: 



X(i) 



Sin7r7 

TT 



tj VE(ii)dii 



or 



P(t) ^ !!!iZ[2 / tr^aE{ti)dti. 

TT Jo 



(24) 



(25) 



In the limit t ^ oo these expressions for X(t) or P(t) are 
MeUin-transforms of the external field, i.e.: 



„ sm 7r7 . ,r„ 
Poo=^ -M[E;j]. 



(26) 



where A^[/;s] = J^ /(t)i''^^di. Thus, the response 
of the CTRW-system to a time-dependent field decays, 
and its polarization tends to a constant value. Interest- 
ingly, the CTRW-system not only ages, but shows a kind 
of "Freudistic" response: the polarization at time t is 
mainly due to the early history of the system, immedi- 
ately after it was prepared in its initial state. 

Let us compare the linear response (polarization vs. 
external field) of a CTRW-system and of a system de- 
scribed by the FFPE. The current through a system de- 
scribed by the FFPE is given by a fractional derivative, 

J(t) - ^^^j\t-hy-'E{t,)dt, (27) 

and the polarization of the system by a fractional integral 

p(i) = i^ / (^ - t2y''nt2)dh. (28) 

For the constant field E(ti) = Eq both expressions, 
Eqs. ( p5|) and (p8|), describe essentially the same time- 
evolution (if one sets (crsin7r7)/7r — a* /T^-j) = ctq), since 
Eqs. ( p5| ) and ( p8| ) are equivalent to each other under the 
change of variable t2 — t — ti. This equivalence doesn't 
hold anymore for the time-dependent-field. Note that 
the main difference of the response, Eq.(28) as compared 
with Eq.(|2F 



is the fact that here the polarization is af- 
fected mainly by the latest events, and that the memory 
of the early history of the system fades away. 



IV. EXAMPLES 

As examples we discuss two simple situations of the re- 
sponse of CTRW and of FFPE systems to external fields. 
We first choose a rectangular pulse switched on at f = t^ 
and off at i = t^ and then a sinusoidal field switched on 
at t = 0. We consider here a highly symmetric system de- 
scribed by a scalar conductivity a^ in which the current 
flows in the direction of the external field. 

For CTRW Eq.(|3|) gives 



J{t) 








L <-. ill! 

t > t^ 



(29) 



which is a causal response concentrated on the time in- 
terval in which the field acts. On the other hand, from 
Eq.(|^ it follows for FFPE that 










(t-t,) 



7-1 



t>t. 



Equations ( |29| ) and (^0|) coincide only if one takes 
tu; = 0, tz —f 00. This limit parallels the findings of 
Ref. [0. In general, however, they are different. Eq.(|29|) 
describes a response concentrated on the time-interval 
of the field-action, t^j < t < tz- no afteraction effects 
are seen. The current never changes sign and has finite 
jumps at i = i„ and t = tz- On the other hand, Eq.(pG) 
shows considerable afteraction: the current does not van- 
ish for t > tz- The current diverges at f = t^ and t — tz 
and changes its sign from positive to negative at t = tz- 
Furthermore, the overall response described by Eq.(pfl) 
is invariant under time translation, i.e. it depends only 
on the differences t — t^ and t ~tz, which is not the case 
for CTRW, Eq.(||). 
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FIG. 1. Shown is the current J(f) in response to two 
rectangular field pulses of unit amplitude acting during the 
intervals 1 < i < 2 and 4 < i < 5. The thick solid line 
reprersents the CTRW-response, while the dashed line repre- 
sents the response of a system described FFPE, see text for 
details. 



Fig.l shows the systems' response to the two pulses of 
unit length and unit amplitude following at times ii = 1 
and ^2 = 4 as follows from Eq.(||) and Eq.(|27|). The 
parameters are cto = 1 and 7 = 1/2. The response to 
the first pulse follows exactly Eqs.(|29|) and (30); hence 
the differences between CTRW and FFPE can be seen 
clearly. The comparison of the response of the two sys- 
tems to a second pulse is also very instructive. For sys- 
tems described by the FFPE the response to the first 
and to the second pulses do not differ (note that the su- 
perposition principle holds and that the responses are 
additive). On the other hand, in the case of a CTRW- 
system the second pulse causes a much weaker reaction 
than the first one. The susceptibility of the system (the 
conductivity) decays with time. 



(30) 



Let us now turn to the response of the system to a sinu- 
soidal force switched on at i = 0, E(t) = £^osin(a;i)0(i). 



The current through the quasiequihbrium (FFPE) sys- 
tem is described by the corresponding fractional deriva- 
tive, J{t) ex qD\^'' E(t). For example, for 7 — 1/2, the 
response to the sinusoidal field is given by: 



^1/2 
J{t) ex -j^Eit) ex cfqEqUJ^/^ [sm{ujt + 7r/4) 

- \/2Gres(v^)j , 



(31) 



whereas 



P{t)(X 



E{t) (X ctqE'ow^^/^ [sin(wt - 7r/4) 






(32) 



In Eqs.(pT|) and (32) Gres(a;) and Fres(a;) are the auxil- 
iary Fresnel integrals; see Ref. |22]. Since both functions 
vanish for x large, the response of the system (J or P) 
at long times tends to be sinusoidal, with a phase shift 
of 7r/4 or — 7r/4 with respect to the acting force. The 
corresponding functions are shown in the upper panel of 
Fig. 2 for (To = -Eq = (^ = 1. 




FIG. 2. The current J(t) (dashed lines) and the polar- 
ization P{t) (full lines) as a response to the sinusoidal field 
E{t) = sint for a system described by FFPE (upper panel) 
and by the CTRW-model (lower panel), see text for details. 

For a CTRW system under the same sinusoidal exter- 
nal field the current through the system decays as 

J{t) ex aQt^'^/^E{t) ex aoEot^'^^^ sincot, (33) 

while the polarization follows: 

P(i) = / J{t)dt(xaoEQ / t-^^^smujtdt (34) 



where S{x) is the Fresnel integral. Note that since 
Ihn^^oD S{x) = 1/2, as time grows the polarization tends 
to a constant value Poo oc Ep/^/uJ as time grows. The be- 
havior of J{t) and P{t), Eqs.(g3|) and (|34| ) is shown in the 
lower panel of Fig. 2 for the same values of parameters as 
before. 



V. CONCLUSIONS 



We have considered the linear response of systems 
governed by CTRW and by FFPE dynamics to time- 
dependent external fields. The form of the response for 
cases described by CTRW displays a "Freudistic" mem- 
ory and no afteraction after switching off the field. This 
differs considerably from the responce shown by systems 
obeying FFPE. 
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